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Abstract 

Background: Gastropod mollusks are one of the most successful animals that have diversified in the fully terrestrial 
habitat. They have evolved terrestrial taxa in more than nine lineages, most of which originated during the 
Paleozoic or Mesozoic. The rissooidean gastropod family Pomatiopsidae is one of the few groups that have 
evolved fully terrestrial taxa during the late Cenozoic. The pomatiopsine diversity is particularly high in the 
Japanese Archipelago and the terrestrial taxa occur only in this region. In this study, we conducted thorough 
samplings of Japanese pomatiopsid species and performed molecular phylogenetic analyses to explore the 
patterns of diversification and terrestrial invasion. 

Results: Molecular phylogenetic analyses revealed that Japanese Pomatiopsinae derived from multiple colonization 
of the Eurasian Continent and that subsequent habitat shifts from aquatic to terrestrial life occurred at least twice 
within two Japanese endemic lineages. Each lineage comprises amphibious and terrestrial species, both of which 
are confined to the mountains in heavy-snow regions facing the Japan Sea. The estimated divergence time 
suggested that diversification of these terrestrial lineages started in the Late Miocene, when active orogenesis of 
the Japanese landmass and establishment of snowy conditions began. 

Conclusions: The terrestrial invasion of Japanese Pomatiopsinae occurred at least twice beside the mountain 
streamlets of heavy-snow regions, which is considered the first case of this event in the area. Because snow 
coverage maintains stable temperatures and high humidity on the ground surface, heavy-snow conditions may 
have paved the way for these organisms from freshwater to land via mountain streamlets by preventing winter 
desiccation in mountain valleys. The fact that the terrestrialization of Pomatiopsidae occurred only in year-round 
wet environments, but not in seasonally dried regions, provides new insight into ancient molluscan 
terrestrialization. 



Background 

Terrestrial invasion is one of the most important events 
in the history of life [1,2]. Terrestrial lineages evolved in 
seven animal phyla, among which vertebrates, arthro- 
pods, and mollusks are representatives of more success- 
ful groups of the epifaunal terrestrial life [3]. Most 
terrestrial lineages are considered to have originated 
under mild or tropical climate during the Paleozoic and 
Mesozoic, whereas few animals became fully terrestrial 
during the Cenozoic [3,4]. In regard to this point, mol- 
lusks are especially unique due to several fully terrestrial 
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and epifaunal lineages that evolved during the Cenozoic 
[5,6], 

Terrestrial invasion of gastropod mollusks has 
occurred in Neritopsina, Cyclophoroidea, Littorinoidea, 
Rissooidea, Ellobioidea, Onchidioidea, Rathouisioidea, 
Succineoidea, and Stylommatophora [7,8], and in parti- 
cular, each of Neritopsina, Rissooidea and Ellobioidea 
has likely achieved land invasion more than once [6,8,9]. 
Most terrestrialization events have occurred during the 
Paleozoic or Mesozoic [2,8-10], and only some members 
of rissooidean families Truncatellidae, Assimineidae, and 
Pomatiopsidae are considered to have colonized to land 
during the Cenozoic [5,6,11,12]. Most truncatellid and 
assimineid snails amphibiously live in intertidal and 
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supratidal zones from brackish water to pelagic areas. 
Terrestrial lineages likely evolved from such ancestors 
[4,6]. In contrast, most pomatiopsid snails live in fresh- 
water habitats [11,13]. Therefore, the process of terres- 
trial invasion in Pomatiopsidae might differ from those 
of former families. 

Pomatiopsidae, which are well known as intermediate 
hosts of Asian schistosomes [14-18], comprises two sub- 
families, Triculinae and Pomatiopsinae. The former 
radiated as aquatic snails in freshwater habitats in 
Southeast Asia [11,19], whereas the latter are distributed 
worldwide (Figure 1) and possesses various lifestyles 
from aquatic to amphibious, littoral, halophilic, and 
even terrestrial [11,20]. The generic diversity of Poma- 
tiopsinae is particularly high in the Japanese Archipe- 
lago, where four of the eight genera including two 
endemics, are recorded (Figures 1 &2) [11]. The diver- 
sity of Japanese Pomatiopsinae is based on two charac- 
teristics, i.e. unique habitats of the endemic genera, and 
their distribution. 

One of the two Japanese endemic genera, Fukuia, live 
amphibiously around rocky walls of steep valleys covered 
with ferns and bryophytes, and moistened by dripping 
water [21,22]. They inhabit only along the mountain 
streamlets, where such habitats are typically found, and 
often occur with pleurocerid freshwater snails [21]. In con- 
trast, the other endemic genus Blanfordia is the sole terres- 
trial group within Pomatiopsidae (Figure 1) [11]. They are 



found on forest floor and tree trunks in coastal to mountai- 
nous areas (Table 1) [11,22], together with fully terrestrial 
gastropod snails such as diplommatinids, camaenids and 
bradybaenids. The distribution of the two genera are lar- 
gely confined to the heavy-snow region facing the Japan 
Sea (Figures 3A-C) [22]. Generally, East Asia is dominated 
by a monsoon climate, characterized by alternating hot 
rainy summers and cold dry winters [23]. The snow- rich 
area is unique in that winter monsoons blowing southward 
from the Japan Sea and high mountain ranges shaping the 
Japanese Archipelago bring enormous snowfalls in winter 
[24,25]. Such a varied climate harbors various unique flora 
and fauna [26,27], and can be closely related to the diversi- 
fication of Pomatiopsinae in the archipelago. 

In this study, we provide the first molecular phylogeny 
encompassing all pomatiopsid taxa occurring on the 
Japanese Archipelago. We also analyzed the relationship 
with continental taxa. In addition, we estimated the tim- 
ing of diversification and terrestrial invasion, which can 
be interpreted using the geological history of Japan. 
Based on these results, we discuss the reason why the 
terrestrial invasion of Pomatiopsidae occurred only in 
the Japanese Archipelago. 

Results 

Phylogenetic analyses 

To investigate the evolutionary history of pomatiopsid 
snails in the Japanese Archipelago, we sequenced the 



Blanfordia (4 s PP .) ^""V- Cecina < 1 s p-) ■ 



Triculinae 

(>20 genera) \ 



Fukuia (3 spp.) 
Pomatiopsis (4 spp.)' 
Oncomelania (2 spp.) 



Idiopyrgus (1 sp.) 



Tomichia (11 spp.) 



d> Terrestrial 
O Littoral 
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Figure 1 A map showing worldwide distribution of Pomatiopsine genera and Triculinae (after Davis [1 1]). Terrestrial taxa occur only on 
the Japanese Archipelago located in East Asia {Blanfordia). Tomichia and Coxiella include several halophilic species occuring on saline lakes. 
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Figure 2 Japanese pomatiopsine snails and their habitats. A. Dripping rocks beside the waterfall covered with bryophytes, inhabited by 
amphibious species F kurodoi. B. Habitats of amphibious species F. multistrioto (beside the waterfall) and terrestrial species B. simplex (forest 
floor). C. Terrestrial species B. joponico on the forest floor. D. Aquatic species F. ooyogii clinging to wet rocky wall. E. Aquatic species 0. minima 
inhabiting mountain streamlet. F. Rice field (dry season) inhabited by seasonally amphibious species 0. h. nosophoro. G. Decaying seaweed 
stranded on the beach, inhabited by C. manchurica. H. F kurodoi niigataensis. I. F. k. kurodoi mating outside the water film. J. F. multistrioto. K. F. 
multistrioto resting on a tree branch. L. B. Integra (terrestrial). M. F ooyogii. N. B. sp. (terrestrial). O. B. bensoni (terrestrial). P. B. simplex mating on 
the forest floor. Q. 0. minimo. R. 0. h. nosophoro. S. C. monchurico found among decaying seaweeds. 



mitochondrial cytochrome oxidase subunit 1 (COl), 16S 
ribosomal RNA genes, and nuclear 18S and 28S rRNA 
regions (Table 1). Maximum parsimony (MP) analysis of 
the combined data resulted in a single most parsimo- 
nious tree (3375 steps, consistency index excluding 
uninformative characters [CI] = 0.3861, retention index 
[RI] = 0.6567). The results of maximum likelihood (ML) 
and Bayesian analyses were largely consistent with the 
relatively well supported clades. Only posterior 



probabilities of P > 0.95 and bootstrap support values of 
P > 75 were considered statistically significant. Hence, 
support values below this significance level are not 
discussed. 

The phylogenetic analyses recovered Pomatiopsidae as 
a monophyletic clade (Figure 4). Triculinae was also 
confirmed as monophyletic. On the other hand, rela- 
tionships among most of the pomatiopsine genera are 
not well resolved, although the pomatiopsine genus 
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Table 1 Sampling information of specimens used in the present study 



Taxon 



Habitat 



Locality 



18S 



28S 



16S 



COI 



POMATIOPSIDAE 
Pomatiopsinae 

Blonfordia bensoni (Adams, 1861) 



Blonfordio Integra Pilsbry, 1924 



Blonfordia japonica (Adams, 1861) 
Blanfordia simplex Pilsbry, 1902 



Blonfordia sp. 

Cecina manchurica Adams, 1861 

Fukuio kurodai kurodai Abbott & 
Hunter, 1949 



Inland forests; terrestrial 



Inland forests; terrestrial, often 
arboreal 



Coastal to inland forests; terrestrial 

Coastal dunes and littoral forests; 
terrestrial 

Inland forests; terrestrial 
Littoral 



Mountain streamlets; amphibious, 
often arboreal 



Fukuio kurodai niigataensis Minato, 
1973 



Fukuio multistriata Abbott & Hunter, 
1949 

Fukuio ooyagii Minato, 1982 

Oncomelania hupensis nosophora 
(Robson, 1915) 

Oncomelania hupensis quadrasi 
(Mollendorff, 1895) 

Oncomelania hupensis robertsoni 
Bartsch, 1946 

Oncomelania minima (Bartsch, 1 936) 



Pomatiopsis lapidaria (Say, 1817) 
Triculinae 

Gammatricula fujianensis (Liu, Zhang 
& Wang, 1983) 

Gammatricula shini (Habe, 1961) 

Gammatricula songi Davis, Chen & 
Yu, 1994 

Lacunopsis sp. 

Neotricula aperta (Temcharoen, 1971) 

Neotricula burchi (Davis, 1 968) 

Roberts iel la sp. 

Tricula bollingi Davis, 1 968 

Triculo fujianensis (Liu, et al., 1983) 

Tricula hongshanensis Tang et al., 
1986 



Mountain streamlets; amphibious, 
often arboreal 



Mountain streamlets; amphibious, 
often arboreal 

Mountain streamlets; aquatic 

Lakes/marshy ground; seasonally 
amphibious 

Lakes/marshy ground; seasonally 
amphibious 

Lakes/marshy ground; seasonally 
amphibious 

Mountain streamlets; aquatic 
Marshy ground; amphibious 

Aquatic 

Aquatic 
Aquatic 

Aquatic 
Aquatic 
Aquatic 
Aquatic 
Aquatic 
Aquatic 
Aquatic 



Nigorikawa, Mori, 
Hokkaido 

Futoro, Setana, 
Hokkaido 

Mt. Mikuni, Shizuoka 

Nou, Itoigawa, Niigata 

Sado, Niigata 1 

Yunohama, Tsuruoka, 
Yamagata 1 

Ayukawa, Fukui, Fukui 

Samani, Hokkaido 

Nanao, Ishikawa 

Betsukai, Hokkaido 

Ohyana, Yurihonjo, 
Akita 

Atsumi, Tsuruoka, 
Yamagata 

Kitamata, Sakata, 
Yamagata 

Shimouchinami, Ohno, 
Fukui 

Takeda, Maruoka, Sakai, 
Fukui 1 

Oyachi, Itoigawa, 
Niigata 

Ikazuchi, Murakami, 
Niigata 

Umeura, Echizen, Fukui 

Iwaya, Mutsu, Aomori 1 
Nirasaki, Yamanashi 

Genbank 

Genbank 

Wajima, Ishikawa 
Sado, Niigata 
Genbank 

Genbank 

Yonaguni Isl., Okinawa 
Genbank 

Genbank 
Genbank 
Genbank 
Genbank 
Genbank 
Genbank 
Genbank 



AB611708* 

AB611712* 

AB611716* 

AB611720* 
AB611724* 
AB611728* 

AB611732* 
AB611736* 
AB611740* 
AB611744* 
AB611748* 

AB611752* 

AB611756* 

AB611760* 

AB611764* 

AB611768* 

AB611772* 

AB611776* 

AB611780* 
AB611784* 



AB611709* 

AB611713* 

AB611717* 

AB611721* 
AB611725* 
AB611729* 

AB611733* 
AB611737* 
AB611741* 
AB611745* 
AB611749* 

AB611753* 

AB611757* 

AB611761* 

AB611765* 

AB611769* 

AB611773* 

AB611777* 

AB611781* 
AB611785* 



AF212906 AY207042 

AB611788* AB611789* 

AB611792* AB611793* 

AF367666 - 

AF212909 - 



AB611710* 

AB611714* 

AB611718* 

AB611722* 
AB611726* 
AB611730* 

AB611734* 
AB611738* 
AB611742* 
AB611746* 
AB611750* 

AB611754* 

AB611758* 

AB611762* 

AB611766* 

AB611770* 

AB611774* 

AB611778* 

AB611782* 
AB611786* 

DQ2 12862 

AF212893 

AB611790* 
AB611794* 
AY676118 



AB611711* 

AB611715* 

AB611719* 

AB611723* 
AB611727* 
AB611731* 

AB611735* 
AB611739* 
AB611743* 
AB611747* 
AB611751* 

AB611755* 

AB611759* 

AB611763* 

AB611767* 

AB611771* 

AB611775* 

AB611779* 

AB611783* 
AB611787* 

DQ1 12287 

AF213339 

AB611791* 
AB611795* 
AF367636 



AF2 12896 AF2 13342 

AB611796* AB611797* AB611798* AB611799* 

EF394879 EF394867 EF394902 

AF212910 - AF212897 AF213343 

AF531540 AY207034 AF531556 AF531541 

AF531543 AY207035 AF531542 AF531544 

AF531549 - AF531548 AF531550 

AF531552 AY207039 AF531551 AF531553 

EF394885 EF394873 EF394893 

EF394888 EF394876 EF394896 
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Table 1 Sampling information of specimens used in the present study (Continued) 



Triculo hortensis Attwood et al., 2003 Aquatic 


Genbank 


- 




EF394883 


EF394871 


EF394900 


Tricula hsiongi Kang, 1984 Aquatic 


Genbank 






EF394889 


EF394877 


EF394897 


Triculo pingi Kang, 1984 Aquatic 


Genbank 






EF394881 


EF394869 


EF394901 


Triculo wumingensis Hu et al., 1994 Aquatic 


Genbank 






EF394884 


EF394872 


EF394892 


ASS I M IN El DAE 














Angustossimineo sotumono Habe, 
1942 


Kikai Isl., Kagoshima 


AB61 


1800* 


AB611801* 


AB611802* 


AB611803* 


Assimineo hirodoensis Habe, 1 942 


Urakami, Nagasaki, 
Nagasaki 


AB61 


1804* 


AB611805* 


AB611806* 


AB611807* 


Poludinellossimineo joponico (Pilsbry, 
1901) 


Hiburi Isl., Uwajima, 
Ehime 


AB61 


1808* 


AB6 11809* 


AB611810* 


AB611811* 


Pseudompholo miyozokii (Habe, 1943) 


Ashikari, Ogi, Saga 


AB61 


1812* 


AB611813* 


AB611814* 


AB611815* 


TRUNCATELLIDAE 














Truncotello pfeifferi Martens, 1860 


Nanao, Ishikawa 


AB61 


1816* 


AB611817* 


AB611818* 


AB611819* 


AMNICOLIDAE 














Akiyoshio koboyoshii Kuroda & Habe, 
1958 


Taga, Shiga 


AB61 


1820* 


AB611821* 


AB611822* 


AB611823* 


LITTORINIDAE 














Echininus cumingi spinulosus (Philippi, 
1846) 


Kikai Isl., Kagoshima 


AB61 


1824* 


AB611825* 


AB6 11826* 


AB611827* 


Littorino pollescens (Philippi, 1846) 


Nago, Okinawa 


AB61 


1828* 


AB611829* 


AB611830* 


AB611831* 



^ype locality 

Taxon names, habitats, localities, and accession numbers are provided. Sequences obtained for this study are marked with an asterisk. 



Oncomelania was sister to Triculinae. Following analysis 
using an approximately unbiased(AU) test [28], the 
monophyly of all pomatiopsine taxa was rejected (P = 
0.049), suggesting that Pomatiopsinae is paraphyletic to 
Triculinae. 

Japanese pomatiopsine taxa were grouped into four 
lineages. Oncomelania and Cecina were both recovered 
as monophyletic, but two endemic genera, Fukuia and 
Blanfordia formed two clades, each of which included 
species from both genera (clades I and II; Figure 4). 
Clade I consisted of F. ooyagii and Blanfordia spp. 
except for B. integra, where F. ooyagii was sister to the 
Blanfordia species. Clade II was composed of B. integra, 
F. kurodai, and F. multistriata. These results indicate 
independent evolution of terrestriality in the two clades. 
Therefore, we tested the alternative hypothesis of the 
occurrence of terrestrial species. The monophyly of all 
Blanfordia species was rejected by approximately 
unbiased (AU) test (P = 0.000001), suggesting that land 
invasion occurred independently at least twice. In clade 
II, B. integra was placed between F. multistriata and F. 
kurodai. However, the detailed branching pattern and 
number of terrestrial invasions were uncertain within 
this clade because the monophyly of B. integra was not 
rejected (AU test, P = 0.07). 

Estimation of divergence time 

To determine the timing of the origin of terrestriality in 
Pomatiopsidae, we estimated the divergence time based 



on the combined phylogeny of four genes. Fossil records 
and previous estimations were used to estimate the 
divergence time of pomatiopsid taxa. For the fossil 
records, Assiminea spp. from the Lower Miocene (23.0- 
16.0 Ma) [29,30] were used for the initial divergence of 
Assimineidae. The root node of "Truncatelloidea" 
(Truncatellidae + Pomatiopsidae + Assimineidae; sensu 
Ponder [31]) was modeled by a Gamma prior (shape = 1 
and scale = 1.5) with values younger than 61.7 Ma 
assigned zero probability, which is justified by Trunca- 
tella minor, the oldest known fossil for this group docu- 
mented from the Lower Paleocene (65.5-61.7 Ma) of 
Belgium [32]. Based on these constraints, a basal split 
within Triculinae was inferred at 16.8 Ma with a 95% 
confidence interval of 21.3-12.2 Ma (Figure 5A). The 
emergence of Japanese aquatic species O. minima and F. 
ooyagii were estimated at 16.0 (21.2-10.6) Ma and 17.4 
(23.7-11.1) Ma, respectively. The terrestrial lineage in 
clade I began to diverge at 6.4 (9.2-3.7) Ma. The diversi- 
fication of clade II began at 7.2 (10.3-4.4) Ma and the 
terrestrial species B. integra occurred during the 
Pliocene. 

We also estimated the divergence time by using the 
split between Oncomelania hupensis quadrasi and O. h. 
robertsoni, which was previously estimated at 3 Ma 
based on allozyme differences [33]. Node ages calculated 
from this estimation were two to three times younger 
than former results. Thus, data from the two calibra- 
tions were not compatible at higher nodes (data not 
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A Cladel 

o • F. ooyagii - Aquatic 
■ B. japonica 

♦ OB. bensoni 
▼ VS. simplex 

* B. sp. 



Terrestrial 




B Clade II (amphibious) 
♦ OF. kurodai 
■ □ F. multistriata 




C Clade II (terrestrial) 
• oS. integra 




D Other Pomatiopsinae 
▼ 0. minima 
a a 0. hupensis nosophora 
+ tyC. manchurica 



130°E 




142°E 36°N 



Figure 3 Maps of Japan showing distributions of Pomatiopsine taxa. (a) Clade 1; (b) amphibious species of clade 2, Fukuio kurodai and F. 
multistriata; (c) terrestrial species of clade 2, B. integra; (d) other pomatiopsine species. Solid and open symbols represent the collecting sites in 
this study and distribution record from the literature, respectively [22,35,72-76]. Solid and dashed curves indicate areas that have more than 100 
mm and 150 mm, respectively, of total precipitation in January (average of 1971-2000 [77]). 
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— G. fujianensis 
G. shini Yonaguni Is." 
G. songi 

Tricula wumingensis 
Neotricula aperta 
Robertsiella sp. 
N. burchi 
fujianensis 



1 00/73/1.01 T. hsiangi 
T. pingi 
T. bollingi 

Lacunopsis sp. 
T. hortensis 
T. hongshanensis 

O. h. quadrasi Philippines 
O. h. nosophora* 
O. h. robertsoni China 
minima Noto* 
O. minima Sado* 
C. manchurica Hokkaido* 
C. manchurica Ishikawa* 
100 /99/1 .or 5 - simplex Fukui* 

B. simplex Yamagata* 
- B. japonica* 
B. bensoni Futoro* 
B. bensoni Mori* 
-S. sp. Samani* 
F. ooyagii 

F. kurodai kurodai Tsuruoka* 
F. k. niigataensis Murakami* 
F. k. kurodai Sakata* 
F. k. kurodai Akita* 
F. k. kurodai Maruoka* 
F. k. niigataensis Itoigawa* 
F. k. kurodai Ono* 
B. integra Itoigawa* 
B. integra Shizuoka* 
F. multistriata* 
P. lapidaria N America 

Angustassiminea satumana 
Pseudomphala miyazakii 



0 



□ J 



Clade I 



Assiminea hiradoensis 

Paludinellassiminea japonica — 
Truncatella pfeifferi — T ru n ca te 1 1 i d a e 



Clade II 

[] Aquatic 
| Amphibious 
| Terrestrial 
0 Littoral 

Assimineidae 



Akiyoshia kobayashii 



Echininus cumingi spinulosus 

Littorina pallescens 



0.05 substitution / site 



Figure 4 A maximum likelihood tree based on the combined dataset of 18S, 28S, 16S, and COI genes (-In likelihood = 20667.0486). 

Numbers on branches indicate maximum parsimony and maximum likelihood bootstrap values followed by Bayesian posterior probabilities. Bars 
on the right of the tree indicate habitat of each taxon. Japanese pomatiopsid taxa are indicated by asterisks. Clades I and II are endemic to 
Japanese Archipelago. 



shown). Although we could not determine which esti- 
mation is more accurate, we employed the former 
results because the node ages in the latter estimation 
were too young to explain the existence of fossils at ear- 
lier periods. 

Climatic factors contributing to distribution of endemic 
genera 

The two Japanese endemic lineages, clades I and II, are 
largely confined to the Japan Sea side of the Japanese 



Archipelago (Figure 3). To determine whether this pat- 
tern is associated with any climate factors, we performed 
logistic regression analysis based on four variables, i.e., 
annual highest and lowest temperatures, and total preci- 
pitations in summer and winter. The logistic regression 
analysis showed significant correlation between winter 
precipitation and the occurrence of the snails, whereas 
contributions of the other variables were not significant 
(Table 2). Although we did not include the amount of 
snowfall due to the significant correlation with winter 
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precipitation (Spearman's rank correlation coefficient, p 
= 0.0054), its difference was also significant between 
inside and outside of the distributions (Wilcoxon's 
signed-rank test, z = -3.009, p = 0.0026). 

Discussion 

Phylogenetic and taxonomic implications 

The phylogenetic analyses supported the monophyly of 
Triculinae. On the other hand, Pomatiopsinae was para- 
phyletic to the former subfamily, and an Asian genus 
Oncomelania became sister to Triculinae (Figure 4). 
These results suggest that taxonomy and phylogeo- 
graphic scenarios should be reconsidered. Traditionally, 
Pomatiopsidae is considered to have diverged into two 
subfamilies in Gondwanaland during the Mesozoic, 
mainly because three pomatiopsine genera shows Gond- 
wanan distribution (Coxiella, Idiopyrgus and Tomichia; 
Figure 1) [11]. However, the close relationship between 
Triculinae and Oncomelania and the estimated diver- 
gence time suggest that Triculinae have evolved more 
recently on the Eurasian Continent, not in Gondwana- 
land. Although we could not provide alternative scenar- 
ios from current phylogenetic results, further study, 
including remaining pomatiopsine genera, is needed to 
reconsider subfamilial taxonomy and worldwide phylo- 
geography of this family. 

Incongruence between current taxonomy and molecu- 
lar phylogeny was also observed on Japanese endemic 
genera Fukuia and Blanfordia. The terrestrial species B. 
Integra was grouped with the amphibious species F. kur- 
odai and F. multistriata, whereas aquatic species F, 
ooyagii was separated from the other Fukuia species 
(Figure 4). Our results indicated the inadequacy of the 
current taxonomy of the two genera, which have been 
primarily based on their habitat [34]. In contrast, inter- 
nal morphologies likely reflect the phylogenetic relation- 
ship. Davis [11] noted that genital structures of B. 
integra are similar to those of Fukuia. Although he 
doubted if Fukuia and Blanfordia are morphologically 
distinct based on such evidence, it is an apparent confu- 
sion resulting from the inadequate generic assignment 
of B. integra, for genital and nervous structures clearly 
differ between species in clade I and II [11,21]. Snails in 
clade II also commonly climb up grasses and shrubs 
during summer (Figure 2K) [11, 35, Y. Kameda, pers. 
obs.]. This behavior characterizes clade II, since tree- 
climbing behavior is known only in these snails and a 
few arboreal species of the assimineid genus Omphalo- 
tropis throughout the superfamily [36]. Hence, B. integra 
is undoubtedly a member of the genus Fukuia. At the 
same time, F, ooyagii should be separated from Fukuia, 
and its generic assignment should be determined 
coupled with the investigation of its soft-part 
morphology. 



Divergence time and patterns of Japanese Pomatiopsinae 

The number of inhabiting pomatiopsine genera is parti- 
cularly high in the Japanese Archipelago (Figure 1). 
Davis [11] considered that the diversity of Japanese 
Pomatiopsinae results from the single colonization and 
following radiation of ancestral Oncomelania in the 
archipelago, and that such rapid evolution presumably 
associated with extensive tectonic activity in the Ter- 
tiary. However, our phylogenetic analysis grouped Japa- 
nese pomatiopsine snails into four lineages, 
Oncomelania, Cecina, and clades I and II (Figure 4). 
These results contradict the scenario of Davis [11], sug- 
gesting that current diversity is partially derived from 
the multiple colonization of pomatiopsine ancestors into 
the Japanese Archipelago. Subsequent habitat shifts may 
have contributed to increasing diversity and produced 
the terrestrial lineages. 

The widespread genus Oncomelania comprises two 
Japanese taxa, O. minima and O. h. nosophora (Figures 
2Q &2R). These two taxa are not closely related, sug- 
gesting that colonization to the Japanese Archipelago 
occurred twice in this genus. The age of divergence 
between O. minima and O. hupensis was estimated 
around 16.0 Ma in the Middle Miocene (Figure 5A), 
suggesting that this divergence was likely invoked by the 
expansion of the Japan Sea, which is thought to have 
occurred during this period [37,38]. O. h. nosophora 
diverged from other continental subspecies more 
recently and likely colonized into Japan during glacial 
periods in the Late Pleistocene [39]. 

The littoral genus Cecina also occurs in the Japanese 
Archipelago. This genus consists of a single widespread 
species, Cecina manchurica, living among the decaying 
seaweed stranded on the seashore (Figures 2G &2S) 
[11]. Although we could not elucidate the divergence 
time of this genus due to low support values, the wide 
distribution was formed by recent dispersions via ocean 
currents based on very low genetic diversity between 
geographically distant populations (Figures 3D &4). 

The divergence between two Japanese endemic 
lineages, clades I and II, was estimated to be around 
21.7 Ma in the Early Miocene (Figure 5A). This period 
corresponds to incipient opening of the Japan Sea, and 
the northern and southern parts of the Japanese Archi- 
pelago that were separated from each other [37,38,40]. 
This geological event might have also split the two 
clades (Figure 5B). The northern part of the Japanese 
Archipelago was further fragmented at later periods 
until the orogenesis of the northeastern Japanese Archi- 
pelago began in the later Late Miocene [37,41]. Specia- 
tion in clade I may be explained by geographic isolation 
due to this geological event, and the narrow distribution 
of each species may reflect the past geographic isolation 
(Figures 3A &5C). Remarkably, B. simplex is widely 
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Figure 5 Maximum clade probability tree displayed as a chronogram and paleogeography around the Japanese Archipelago (A) 

Divergence times calculated in BEAST based on fossil records (outgroups not shown). Japanese pomatiopsid taxa are indicated by asterisks. 
Horizontal bars in the tree give the 95% credible interval for the age of each node. Numbers above and below branches represent mean node 
ages and posterior probabilities, respectively. (B, C) Hypothesized paleodistribution of clades I and II. Paleogeographic maps of East Asia were 
modified from lijima &Tada [37]. Shaded areas indicate land. (B) Early Miocene (23-18 Ma), divergence of two clades. (C) Middle to late Pliocene 
(3.5-2 Ma), ongoing diversification and range expansion within each clade. 
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Table 2 Chi-square statistics from the logistic regression 
analysis 



Variables 


x 2 


P 


Highest temperature 


1.6651 x 10~ 9 


1 


Lowest temperature 


9.7577 x 10~ 7 


0.9992 


Summer precipitation 


2.9000 x 10~ 8 


0.9999 


Winter precipitation 


34.7674 


< 0.0001 



Degrees of freedom was 1 for all variables. 



distributed along the Japan Sea and occurs on coastal 
dunes as well as littoral forests, whereas the other spe- 
cies inhabit rather inland forests [22,42]. The low 
genetic differentiation and habitat of B. simplex suggest 
that a geological barrier has been absent for the coastal 
species. 

The southern part of the Japanese Archipelago, where 
clade II was thought to have inhabited (Figure 5B), 
remained as a relatively large landmass. The land area 
had expanded northward in the Late Miocene along 
with the beginning of orogenesis in the northeastern 
part of the archipelago [37,38]. The diversification of 
clade II likely started in the Late Miocene (around 7.2 
Ma; Figure 5A), indicating that tectonic activities might 
have enabled the range expansion of clade II to north- 
eastern Honshu (Figure 5C). 

Routes of terrestrial invasion 

The terrestrial pomatiopsid species are found in clades I 
and II (Figure 4). The snails in both clades are largely 
restricted to mountains in the Japan Sea side of the 
Japanese Archipelago (Table 1; Figures 3A-C), where 
winter precipitations are especially high among the dis- 
tributions of pomatiopsine genera (Figure 6). Because 
winter precipitation falls as snow in this region, terres- 
trial invasion may have occurred via mountain stream- 
lets that are covered with heavy snowfall in winter. This 
evolutionary hypothesis contrasts with the current 
knowledge of terrestrial invasion in two points. First, 
evolution of terrestriality most likely occurred via moun- 
tain streamlets. Generally, that animals would evolve 
terrestriality near rivers and streams than in lakes and 
swamps seems less likely due to adaptations for air 
breathing, particularly in the latter environment with its 
frequent reductions in dissolved oxygen [3]. Second, ter- 
restrial invasion likely occurred in a cool-temperate 
region. Colonization to land is considered to have 
occurred primarily in the tropics during mild climate 
periods [3,4]. To our knowledge, this is the first case in 
which terrestriality might have evolved in a cool region. 
Therefore, the questions on how terrestriality was 
acquired beside the streamlets and which climatic fac- 
tors enabled this evolution under such cold climates 
remain to be elucidated. 
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Figure 6 Monthly precipitations in the distributions of 
pomatiopsine genera. Values for each genus, except for Idiopyrgus 
and Tomichia, are average of 5 localities which was randomly 
selected from the distribution. The precipitation data were obtained 
from the World Meteorological Organization (WMO; http://www. 
wmo.int/). Blue and orange lines represent genera in Northern and 
Southern Hemisphere, respectively. Distribution of Fukuia and 
Blanfordio has especially high precipitation during winter (shaded 
period). B, Blanfordia; Co, Coxiella; F, Fukuia; I, Idiopyrgus; O, 
Oncomelonia; P, Pomotiopsis; T, Tomichia. 

K J 

The aquatic and amphibious species of clades I and II, 
which likely represent ancestral habitats of terrestrial 
taxa, inhabit mountain streamlets and rocky walls that 
are constantly moistened by dripping water and covered 
with ferns and bryophytes (Figure 2) [21,22]. When 
these habitats become drier, both amphibious and aqua- 
tic species are able to resist the short-term drying by 
crawling into crevices and wet litter [21]. Correlation 
between desiccation tolerance and wetness of habitat is 
observed in some animals [11,43,44]; thus, frequent 
droughts may have promoted the ancestral freshwater 
snails to the more amphibious lifestyle beside the moun- 
tain streamlets. Amphibious snails in some populations 
can actively move onto terrestrial substrata and mate 
outside the water (Figure 21). Therefore, the colonization 
to terrestrial habitats was presumably achieved by 
becoming more tolerant to desiccation. The availability 
of dissolved oxygen might not have influenced the habi- 
tat shifts of these snails because their ancestors may 
have preadapted to amphibious habitats [11,12]. 

The snails in clades I and II mainly occur on the 
Japan Sea side of the Japanese Archipelago (Figure 3), 
whose climate is unique in East Asia. This area is gener- 
ally dominated by a monsoon climate, characterized by 
alternating hot rainy summers and cold dry winters 
[23]. However, the Japan Sea side of the archipelago has 
exceptionally high winter precipitation due to the winter 
monsoon blowing southward from the Japan Sea and 
the high mountain ranges that bring enormous snowfall 
[24,25,45]. This region harbors various unique flora and 
fauna that are known as "Japan Sea elements" [26,27], 
whose distributions clearly correspond to the amount of 
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snowfall [46-48]. These organisms are not always resis- 
tant to the stress of winter coldness and desiccation, 
and often need snow coverage under which stable tem- 
peratures around 0°C and high humidity are maintained 
[49-54]. Snails of clades I and II are also considered as 
Japan Sea elements [22,35]. The logistic regression ana- 
lysis showed significant correlation between winter pre- 
cipitation and occurrence of the snails (Table 2), 
suggesting a high degree of need for high winter precipi- 
tation, which brings snow coverage, to ensure their 
survival. 

Heavy-snow environments in the Japanese Archipelago 
were established by the formation of Japan Sea and high 
mountain ranges in the archipelago [55,56]. The orogen- 
esis of the northern Honshu occurred since the late 
Miocene [37,38,41], whose period corresponds to the 
diversification and terrestrial invasion of Japanese ende- 
mic clades (Figure 5A). According to the vegetation in 
the late Miocene to Pliocene, northern Honshu was 
cool-temperate and later became colder [57,58]. The 
amount of snowfall increased with the growth of high 
mountain ranges during the Pleistocene [56]. Although 
the amount of snowfall during these periods is uncer- 
tain, arid or semi-arid climates did not appear in this 
region [24,56], suggesting the presence of constantly 
moist environments that were suitable for the survival 
of snails in clades I and II. Therefore, the establishment 
of snowy regions may have paved the way for these 
organisms from freshwater to land via mountain stream- 
lets by preventing winter desiccation in mountain 
valleys. 

The fact that terrestrial invasion in Pomatiopsidae did 
not occur except in the snowy region provides an indi- 
cation of molluscan terrestrialization. Many pomatiopsid 
snails lived in the floodplains of continental Asia [11], 
where the dry and wet seasonal cycle is dominant (Fig- 
ure 6). However, such climate seems insufficient to lead 
the seasonal amphibious snails to terrestrial habits. In 
contrast, year-round wet environments were most likely 
necessary for amphibious snails to become terrestrial, 
contrary to the belief that seasonal drought accelerates 
terrestrial invasion [3]. The terrestrial lineages, especially 
in clade II, have specific behaviors when crawling up 
vegetation [11,35]. Ancestral amphibious lineages in 
water films found on wet rocky walls would have 
enabled the preadaptation of this type of behavior 
accompanied by a broad, strong foot. 

Conclusions 

The generic and ecological diversity of Pomatiopsinae is 
particularly high in the Japanese Archipelago, with its ter- 
restrial taxa occurring only in this region. Our phyloge- 
netic analyses suggested that the Japanese Pomatiopsinae 
derived from the multiple colonization from the Eurasian 



Continent and subsequent habitat shift to the terrestrial 
life occurred at least twice within the Japanese endemic 
lineages (clades I and II; Figure 4). The evolution of ter- 
restriality possibly occurred beside mountain streamlets 
of snowy regions, where snow coverage shelters the 
ground surface from winter desiccation and coldness 
[51,54]. The estimated divergence time also suggested 
that diversification of the endemic clades began in the 
Late Miocene, which is the period corresponding to the 
orogenesis of the Japanese Archipelago and expansion of 
the Japan Sea (Figure 5). The emergence of the heavy- 
snow region is associated with these tectonic activities 
[56]; therefore, the terrestrial invasion of pomatiopsid 
snails may have been promoted by the formation of the 
archipelago and the establishment of snowy conditions 
that led the way from freshwater to land via mountain 
streamlets. The fact that terrestrial invasions of Poma- 
tiopsidae occurred only in year-round wet and cool 
environments, but not in seasonally dried regions, gives 
new insight into molluscan terrestrialization. 

Methods 

Sampling 

To investigate the evolutionary history of Japanese 
pomatiopsids, we collected 13 taxa from 23 localities to 
represent all Japanese pomatiopsids. The sequence data 
of non-Japanese pomatiopsids were obtained from Gen- 
Bank, but species that lack two or more examined gene 
regions were excluded from the phylogenetic analysis. 
As a result, 13 triculine species and 3 non-Japanese 
pomatiopsine taxa were included in the analysis, 
although they do not encompass representatives of all 
pomatiopsine genera. However, it would not matter in 
discussing the evolution within Japanese taxa because 
excluded pomatiopsine genera are unlikely related to 
our focusing taxa based on genital structures [11]. For 
outgroups, we used four assimineids, one truncatellid, 
one amnicolid, and two littorinids. Detailed information 
on examined samples is shown in Table 1. 

Molecular methods 

Total DNA was isolated following a previously described 
method [59]. We sequenced the partial region of the 
nuclear 18S rDNA and 28S rDNA, mitochondrial cyto- 
chrome oxidase subunit 1 (COl), and 16S rDNA. Poly- 
merase chain reaction (PCR) conditions and primers 
used are shown in Additional File 1. PCR products were 
purified using Exo-SAP-IT (Amersham Biosciences, Little 
Chalfont, Buckinghamshire, UK). Sequencing was per- 
formed using a BigDye™ Terminator Cycle Sequencing 
Ready Reaction Kit (Applied Biosystems, Foster City, CA) 
and electrophoresed on an ABI 3130 sequencer (Applied 
Biosystems). The obtained sequences were deposited in 
the DNA Data Bank of Japan (DDBJ) (Table 1). 
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Phylogenetic analyses 

The alignment of the COl gene required no gaps; 16S, 
18S, and 28S sequences were aligned using MUSCLE 
v3.8 [60]. To eliminate the uncertainty of the alignments 
in these gene regions, GBLOCKS v0.91b [61] was used 
to select those regions of the aligned sequences that are 
confidently aligned for analysis (Additional File 2). Phy- 
logenetic trees were obtained by maximum parsimony 
(MP), maximum likelihood (ML), and Bayesian methods. 
MP analysis was performed using PAUP* 4.0bl0 [62] 
and heuristic searches were conducted with 10 random 
addition analyses using equal character weights and 
tree-bisection-reconnection (TBR) branch swapping. 
Prior to ML and Bayesian analyses, we used the pro- 
gram Kakusan3 v3. 5. 2009. 11.03 [63] to select appropri- 
ate models of sequence evolution (Additional File 2). 
Based on the selected models, ML analysis was per- 
formed with 1000 iterations of the likelihood ratchet 
[64,65] using TreeFinder [66] and Phylogears version 
1.5.2009.12.29 [67]. Nodal support for the MP and ML 
analyses was assessed using bootstrap analyses with 
1000 replications. Bayesian analysis was performed using 
MrBayes v3.1.2 [68,69]. The analysis was carried out by 
two simultaneous runs, each of which consisted of run- 
ning four simultaneous chains for 10 million generations 
and sampling trees every 1000 generations. We dis- 
carded the first 1001 trees as burn-in and the remaining 
18000 samples were used to estimate the tree topology, 
branch length, and substitution parameters. The align- 
ments and obtained trees were deposited in TreeBASE 
http://www.treebase.org/treebase-web/home.html under 
accession number S11213. 

Approximately unbiased (AU) test 

Given that the phylogenetic analyses indicated that the 
terrestrial genus Blanfordia and aquatic/amphibious 
genus Fukuia formed mixed clades (see Results), we per- 
formed the AU tests [28] to test an alternative hypoth- 
esis for the evolution of terrestrial species. The 
alternative trees were obtained by the likelihood ratchet 
method as described above under the topological con- 
straint that all Blanfordia species were monophyletic (i. 
e., terrestrial invasion occurred once) and that B. Integra 
was monophyletic, indicating that terrestriality evolved 
at least twice. In addition, the AU test was also per- 
formed to test the monophyly of Pomatiopsinae, because 
it was not recovered as monophyletic and because sup- 
port values among pomatiopsine genera were low. The 
AU test was conducted based on 1 million replications 
using TreeFinder. 

Estimation of divergence time 

Approximate divergence times were estimated using a 
relaxed clock method [70] as implemented in BEAST 



1.5.3 [71]. Divergence times were calculated using 
previous estimations and fossil records. We first used 
the split between Oncomelania hupensis quadrasi and 
O. h. robertsoni, which was previously estimated at 3 
Ma based on allozyme differences [33]. For the fossil 
records, three fossils were used to place priors on 
node ages. A normally distributed estimate prior of 
19.5 ± 3.5 Ma was used for the basal divergence of 
Assimineidae, based on Assiminea from the Lower 
Miocene (23.0-16.0 Ma) of Florida, USA [29] and the 
Marshall Islands [30]. The basal splitting of "Trunca- 
telloidea" (Truncatellidae + Assimineidae + Pomatiop- 
sidae; sensu Ponder [31]) was modeled by a Gamma 
prior (shape = 1 and scale = 1.5) with values younger 
than 61.7 Ma assigned zero probability, which is justi- 
fied by Truncatella minor, the oldest known fossil for 
this group documented from the Lower Paleocene 
(65.5-61.7 Ma) of Belgium [32]. Divergence times for 
the remaining nodes in the phylogenetic tree were 
estimated with BEAST software. The Yule process was 
used to model the speciation process. The Monte 
Carlo Markov chain was run eight times for 10 mil- 
lion generations sampled every 1000 generations to 
ensure that the effective sample size (ESS) values were 
above 200 for most of the parameters. The first 1001 
trees of each run were discarded as burn-in and the 
remaining trees were combined to produce an ultra- 
metric consensus tree using LogCombiner and 
TreeAnnotator vl.5.3 (included in the BEAST soft- 
ware package). 

Climatic factors contributing to distribution of endemic 
genera 

The two Japanese endemic genera, Fukuia and Blan- 
fordia, are largely confined to the Japan Sea side of 
the Japanese Archipelago, i.e., heavy-snow regions 
(Figure 3). The distributions of these genera seem to 
be restricted by the amount of snowfall, but other cli- 
matic factors may have affected the survival of these 
snails. Therefore we performed logistic regression ana- 
lysis to investigate the association between climatic 
factors and the occurrence of these snails. We used 
annual highest and lowest temperatures, and total pre- 
cipitations in summer (between June and August) and 
winter (between December and February) as climatic 
factors. The climatic data were obtained from the 
Japan Meteorological Agency http://www.jma.go.jp. 
Within the distribution of Fukuia and Blanfordia, 14 
localities were selected to cover the entire range of 
these genera (Additional file 3). For the localities out- 
side the distributions, we selected 11 weather stations, 
each of which is nearest to an inside-distribution 
locality across the high mountain ranges and has simi- 
lar altitude. 
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Additional material 



Additional file 1: Information on primers and PCR conditions used 
in this study. 

Additional file 2: Information on sequence alignments and models 
of sequence evolution for maximum likelihood analysis. 

Additional file 3: Information on localities used for the investigation 
of climatic factors. 
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